Single parameter scaling in 1-D Anderson localization. Exact analytical solution 
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The variance of the Lyapunov exponent is calculated exactly in the one-dimensional Anderson 
model with random site energies distributed according to the Cauchy distribution. We derive an 
exact analytical criterion for the validity of the single parameter scaling in this model. According to 
this criterion, states with energies from the conduction band of the underlying non-random system 
satisfy single parameter scaling when disorder is small enough. At the same time, single parameter 
scaling is not valid for states close to band boundaries and those outside of the original spectrum 
even in the case of small disorder. The obtained results are applied to the Kronig-Penney model 
with the potential in the form of periodically positioned (5-functions with random strengths. We 
show that an increase in the disorder can restore single parameter scalingbehavior for states from 
band-gaps of this model. 



I. INTRODUCTION 

The hypothesis of single-parameter scaling (SPS) in 
the context of transport properties of disordered conduc- 
tors was introduced in Ref. |l|. It was suggested that 
scaling properties of the conductance g are determined 
by a single parameter, the conductance itself, through a 
scaling equation 



d(lnL) 



Pig), 



(1) 



where L is the size of a sample. The nature of the scal- 
ing conductance, g was debated for some time until it 
was understoodH that scaling in the theory of localiza- 
tion must be interpreted in terms of the entire distribu- 
tion function of conductivity rather then in terms of its 
momentums. SPS in this case means that the distribu- 
tion function of g is fully determined by a single param- 
eter, which obeys scaling equation (|l]). In Ref. which 
was concerned with scaling properties of one-dimensional 
disordered conductors, Anderson, et al. proposed the pa- 
rameter 



(2) 



as a scaling parameter suitable to describe the fluctua- 
tions of the conductivity. In the limit i — s- oo the in- 
troduced parameter takes a non-random value, 7, which 
is the inverse localization length, lioc, or the Lyapunov 
exponent (LE) characterizing the spatial distribution of 
electron's wave functions .Bu It was suggested in Ref. ^ 
that the introduced parameter has a normal distribution 



and does not exhibit anomalously large fluctuations. Cal- 
culations carried out in Ref. ^ showed that the variance 
of LE, cr^, scales according to the law of large numbers 
cr^ ^ 1/L, and is related in a universal way to LE:B 

a' - l/L. (3) 

This relation is the essence of SPS in the case of strong lo- 
calization, as it presumes that two parameters of the nor- 
mal distribution of LE are reduced to one scaling param- 
eter, 7. According to the scaling theory (see Ref. |l] and 
references therein) almost all states in one-dimensional 
systems are localized, and these systems are, therefore, 
always in the regime of strong localization in the asymp- 
totic limit L ^ lioc- 

Eq. was first derived in Ref. || within the approx- 
imation known as the random phase hypothesis, which 
assumes that there exists a microscopic length scale over 
which phases of complex transmission and reflection co- 
efficients become completely randomized. Under similar 
assumptions, Eq. (|^) was rederived lateiJpy several au- 
thors for a number of different models .Bilil 

Landaucr's representation of the conductance in 
terms of tra.ns.m.ission coefficients for different scatter- 
ing channelsE3E3 reduces the study of the conductance 
in quasi-one-dimensional wires to the analysis of scat- 
tering orj-transfer matrices. Within the transfer-matrix 
approachta the problem is further reduced to the study 
of statistical properties of the products of random matri- 
ces (see Ref. |lj and references therein). In this context, 
the self-averaging of LE and its normal distribution in 
the asymptotic limitpL|-^ 00 are rigorously established 
mathematical facts.BEJO SPS expressed by Eq. (0) was 
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also established in the limit of strong localization for a 
quasi-one-dimensional geometry in Ref. |l^,|l9| with .-the 
use of the Dorokhov-Mello-Pereyra-Kumar equationBj 

The SPS hypothesis has also been verified in the regime 
of weak localization, which exists in the conducting phase 
of three dimensional conductors and in the limit of a 
large number of scattering channels in the quasi-one- 
dimensional geometry. In contrast with the case of strong 
localization, the distribution of the conductance (rather 
than the logarithm of the conductance) has the Gaus- 
sian form with two independent parameters. The vari- 
ance of the conductance, however, was found to be a 
universal numbeiOc3 leaving one again with a single 
scaling parameter - the average conductance. The uni- 
versal conductance fluctuations, which were first discov- 
ered in three-dirpepaional conductors within a diagram- 
matic approach ,LirL3 were later reinterpreted from the 
point of view of the random matrix theory in a quasi- 
one-dimensional geometry (see for review Ref. |l5| ). 

The scaling properties of conductivity have also been 
studied numerically by a number of authors. The general 
concepts of the scaling theory expressed by Eq. (|^) were 
verified by means of Green's function and transfer-matrix 
approaches generalized for two- and three-dimensional 
systems in Ref. The log-normal distribution of the 
conductance and SPS Eq. (^ has been confirmed for 
the Anderson model (AM) in numerical simulations in 
Ref. In the one-dimensional situation, the existence 
of SPS has been alspjobtained in simulations of AM with 
correlated disorder £3 and scalar wave propagatipHia su- 
perlattices with different models of randomness E3^c£l 

The zero energy state in one-dimensional models with 
off-diagonal disorder (random hopping models) repre- 
sents a special case. These models demonstrate |a-|delo- 
calization transition in the vicinity of zero energyn3 con- 
trary to the conclusion of the scaling theory that such a 
transition is absent in one-dimensional systems. The SPS 
relation (^ between the standard deviaJ-ion of LE and its 
mean value is also violated in this case.L3 Unusual proper- 
ties of this model are due to a so called chiral symmetry, 
which is characteristic of the state with E — in this 
model. In further discussion we will ignore this special 
case and refer to regular situations, which include models 
with diagonal disorder and random hopping models away 
from the critical E — point. 

Simultaneously with numerous confirmations of the ex- 
istence of SPS, the limits of its validity have been the 
subject of intensive discussions (see, for example, Ref. 
p9| and references therein). As we mentioned above, the 
original condition for SPS, postulated in Ref. |^, invokes 
the hypothesis of the phase randomization. This hypoth- 
esis implies that the phases of complex transmission and 
reflection coefficients become completely randomized at 
the distances much smaller than the localization length. 
Phase randomization was numerically studied for AM in 
Ref. Il^, where it was shown that for small disorder the 
phases indeed become uniformly distributed at a scale 
much shorter than the localization length. This does not 



happen, however, for states in the center of the original 
conduction band. Numerical calculation of AM in Ref. ^ 
and analytical consideration of a model with periodically 
positioned scatterers in Ref. ^ obtained a nonuniform 

thase distribution for such states. It was shown in Ref. 
I that for the states at the center of the band there ex- 
ists a new length, a phase relaxation length, l^. As soon 
as the length of the sample exceeds l^p, the phase dis- 
tribution approaches a stationary but nonuniform form. 
Under certain conditions, the relaxation length rather 
than the mean-free-path was found to determine the lo- 
calization length. At the same time, neither analytical 
nor numerical studies of the states at the center of the 
conduction band found violations of SPS. These results 
cast doubt upon the relevance of phase randomization for 
SPS. In this paper we show that the condition for SPS is 
not phase randomization. 

When local disorder is strong, thp_phase distribution 
was found to never become uniform,Ej and the probabil- 
ity distribution of LE, ip-this case, is controUecLby two in- 
dependent parameters.OLj It was indicatedcj however, 
that even in the case of an extremely nonuniform phase 
distribution, the deviations from SPS are rather limited. 

The hypothesis of phase randomization lies at the foun- 
dation of all existing theoretical approaches to statistical 
properties of conductaace, including those based upon 
random matrix theoryJlj An additional requirement cru- 
cial for Eq. d) can be called "local weakness of dis- 
order. "pJn_^calculations based upon the random matrix 
theory Mr 1 this requirement is set as a limit when the 
cross-section of each individual scatterer tends to zero, 
while the density of the scatterers tends to infinity keep- 
ing the localization length constant. It is commonly be- 
lieved that in the regime of strong localization (lioc *C L) , 
SPS holds provided that the local disorder is weak, so 
that the localization length exceeds all microscopic length 
scales of the model. Increase of the disorder leads to re- 
duction of the localization length, and eventually violates 
SPS. 

Results which apparently contradict this well estab- 
lished understanding of the crossover between SPS and 
statistics with two independent parameters, were recently 
reported in Ref. The system considered in Ref. ^ be- 
longs to the class of Kronig-Penney-like models (KPM), 
which have been intensively studied (see, for example, 
Refs. and references therein). The original spec- 

trum of KPM contains multiple bands separated by band- 
gaps. Disorder not only localizes states within the orig- 
inal pass bands but -also creates tails of localized states 
in former band gaps.cl According to Ref. |3^ the spectrum 
of the system is divided into two groups of states with 
different scaling behavior: SPS holds for states from the 
conduction bands of the initial spectrum, and is violated 
for states from initial band-gaps. Moreover, this violation 
of SPS for band-gap states occurs even for weak disorder, 
and turns out to be much more dramatic, than the phase 
randomization approach would predict £3 

Occurrence of states outside of the initial conduction 
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bands is known to be a model independent phenomenon. 
It seems plausible, therefore, that the coexistence of SPS 
and non-SPS states found in Ref. |3^ is not a particular 
feature of KPM but rather a general property of quantum 
disordered systems. 

The main objective of the present paper is to re- 
examine the problem of scaling properties of conductance 
in one-dimensional systems and to derive SPS Eq. (|^) 
without the assumption of phase randomization. This 
calculation allows us to formulate a "correct" criterion 
for SPS, and to understand the nature of its violation 
reported in Ref. p2[ The main results of this paper were 
outlined in Ref. 

The paper is organized as follows: In Sec. II we formu- 
late the model within which we calculate the variance of 
LE. The details of the calculations are presented in Sec. 
III. The new criterion of SPS is derived and analyzed in 
Sec. IV. In that section we also complement our ana- 
lytical calculations with numerical simulations of a more 
generic model. Comparison with the latter helps us to 
distinguish between universal features of our results and 
those specific to the selected model. The transition be- 
tween SPS and non-SPS states is discussed in Sec. V. 
We conclude in Sec. VI. 



solution, we assume that parameters e„ or Vn are dis- 
tributed with the Cauchy probability density (the Lloyd 
model) : 



Pcix) = 



1 



TT + (.T — Xq)^ ' 



(7) 



where xq — or Vq for AM or KPM, respectively. Pa- 
rameters xq and r represent the mean value of the ran- 
dom variable x and the width of the distribution, respec- 
tively. Although r characterizes the strength of disorder 
in the system, it cannot be interpreted as a second mo- 
ment of the distribution (0), because the latter does not 
exist. AM with the probability distribution ^ Hr^^ '^^ 
the first models where LE was evaluated exactly.cl'c3 The 
probability distribution of parameters J7„, which enter 
equations of motion (^, has the same form as Eq. (|^) 
with the following parameters: 



{Un 



Uo 



2 cos(fca) + y sin(fca), KPM 
E, AM 



(8) 



II. THE DESCRIPTION OF THE MODEL 



|sin(fca)| 



KPM 
AM 



(9) 



Let us consider a one-dimensional tight binding model 
with diagonal disorder, which is described by the follow- 
ing equations of motion 



-0„+l -I- tpn^l - Unl/Jn = 0, 



(4) 



In the absence of disorder, the energy spectrum of the 
model is determined by the condition: \Uo\ < 2. In AM 
this leads to a single conduction band ~2 < E < 2. In 
KPM there exist multiple bands separated by band-gaps. 
Allowed values of the energy variable belong to intervals 



where ■(/;„ represents the wave function of the system at 
the n-th site. In Eq. (jj) the hopping integral is cho- 
sen to be equal to one, so it sets the energy scale in the 
system. The concrete meaning of t/„ depends upon the 
interpretation of the model (m. There are two appar- 
ently different models that can be described by Eq. (||). 
In the first, this equation represent a classical AM, with 
Un defined as 



Un = -E 



(5) 



where E is the energy of a particle, and e„ is the random 
site energy. Second, it can be shown (see, for examplc,Q) 
that the Schrodinger equation for KPM with a random 
potential formed by periodically positioned (5-functions 
with random strengths, Vn, also reduces to the form (^) 
with ip„ being the values of the eigenfunctions at the sites 
occupied by the (5-potentials. In this case, C/„ is defined 
as follows :□ 



Un 



2 cos(A:a) H — sm{ka)^ 



(6) 



where A; = v i5 is the energy variable and a is the period 
of the structure. To be able to obtain an exact analytical 



k„ < k < 7rn, 



where k'^ obeys the equation 



1,2,3.. 



tan 



tan 



Vo 



n odd 



(10) 



(11) 



(12) 



The high energy bouadaries of each band correspond to 
so called resonanceaa because disorder does not effect 
transport at these particular energies. This fact can eas- 
ily be seen from Eq. (^), where Tu for KPM becomes 
zero for all k — im. The presence of these resonances 
is a specific property of the model under consideration 
caused by the strict periodicity in the positions of site 
potentials. Similar resonances arej-also present in other 
models such as t]ae,|-dimer modeCJ or models of ran- 
dom superlattices.E3E3 The resonances dissappear once 
one destroys the exact periodicity in the positions of 5- 
functions or allows for random variations in the width of 
superlattice's layers. 
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The main objects of our study are the finite size LE, 
7(L), and its variance cr^. 7|(|^) can be defined for the 
model under consideration asl 



j{L) = jlnrN, 



(13) 



where N = L/ais the total number of sites in the system 
and r^r is the envelope of the wave function 



rN 



.1/2 



(14) 



As we discussed in the Introduction, "f{L), takes a non- 
random value, 7, in the limit L ~* oo. This limiting value 
can also be considered as aMpaverage of 7(L) over different 
realizations of the system.Qu For large but finite L, 'y{L) 
exhibits finite size fluctuations whose distribution func- 
tion asymptotically approaches the jGLapasian form with 
the variance, cr^, decreasing as l/L.aOO 

The average LE, 7, in the considered model was first 
calculated in Ref. |3^. It turns out that the method devel- 
oped in that paper (see also Ref. ^ can be as well used 
for exact calculation of a^. The method is based upon 
the representation of LE in terms of the phase variable 
Zn, defined as z„ = ipn/^n-i, which obeys the following 
equation of motion 



Zn + Z, 



n-1 



Finite size LE can be expressed in terms of z„: 

N 

L 



7 



1 ^ 

(L) = -^ln|z„| 



1 



2L 



In 



1 



z 



N+l 



0" 



(15) 



(16) 



If Zn is a stationary random function of n, that is, a dis- 
tribution of Zn is independent of n, the first term in Eq. 
( |l6| ) is of the order of unity while the second term is of 
the order of (1/i) and disappears in the limit L — > 00. 
The expression for LE, therefore, takes the following form 



1 ^ 

lim — In |z„| (In |z„|) 

L— »oo Li ^ — ^ 



(17) 



where the average on the right-hand side is taken over 
the stationary distribution of z. 

The asymptotic expression for variance of LE can be 
obtained from Eq. ([iq): 



1 ^ 

o-^ = — ^ [(lnz„ilnz„) - (Inz™) (lnz„)] 



(18) 



m,n— 1 



and is valid as long as the system's size L is much greater 
than the correlation radius of z„, which we assume to be 
finite. 



III. VARIANCE OF THE LYAPUNOV 
EXPONENT IN THE LLOYD MODEL 

A. Two-point distribution of the phases Zn 

Calculation of the variance from Eq. ( [T8| ) re- 
quires knowledge of the two-point distribution function, 
Piizn, Zm), of the phases z. Our calculations of this func- 
tion are based upon representation of a joint distribution 
of multiple random variables as the product of marginal 
and conditional distributions: 



P2{Zn,Zn+k) = Pl{Zn)P{Zn\Zn+k), 



(19) 



where Pi{zn) is a stationary probability distribution of 
Zn, and P{zn\zn+k) denotes a conditional probability dis- 
tribution of Zn+k provided that z„ is fixed. With the help 



of Eq. (15) the latter probability can be written as 



P{Zn\Zn+k) = / <5 {zn+k + ^^n+fe-i - Un+k-l) 



X P{Zn\Zn+k-l, Un+k-l)dUn+k-ldZn+k-l, 

where P{zn\zn+k-iTUn+k-i) is a joint probability of 
Zn+k~i and Un+k~i- It follows from the structure 
of Eq. (|l^) that z„ depends only upon values of 
the random parameter Um at preceding sites m < n, 
and thus is independent of [/„. The joint probability 
P{zn\zn+k-i,Un+k-i), therefore, can be factorized and 
integration over Un+k-i can be carried out. The result is 
the following recurrent relation between P{zn\zn+k) and 

P{Zn\Zn+k-l)- 



P{Zn\Zn+k) = / P {Zn\Zn+k-l) 



(20) 



Pc{ Zn+k + ^„+fe_i) dZn+k- 



where Pc is the Cauchy distribution introduced in Eq. 
(0). The advantage of the Cauchy distribution is that 
recurrence ( pO| ) can be solved exactly. The conditional 
probability obtained has again the form of the Cauchy 
distribution, which can be conveniently presented in the 
form 



P{,Zn I Zn 



1 



{Zn+k - Cfe) {Zn+k - Ck)' 



(21) 



where the asterisk denotes complex conjugation and pa- 
rameters ^fe obey the following equation 



(22) 



Eq. (|2l|) for P{zn\zn+k) and Eq. (||) for ^j. have ex- 
actly the same form as those obtained in Ref. ^ for 
the one-point distribution Pi(z„). However, in the case 
of the one-point distribution one looks for a stationary 
solution of Eq. (p^), while the conditional distribution 
P{zn\zn+k) rcquircs that Eq. ( |2^ ) is solved with the ini- 
tial condition 



Co 



(23) 



4 



This solution can be presented as 

gk _ g-k _ (^gk+i _ 



-k-l\ 



gk-l _ _ (^fc _ ^-fe) 



(24) 



where S is the fc- independent solution of Eq. (|22|), which 
obeys the stationary version of Eq. (^2|) 



6 + S- 



Uq + ir. 



(25) 



Real and imaginary parts of S determine the center and 
the width of the one-point distribution Pi(z„) 



PliZn) 



luiS 



1 



TT (z„ - (5) (z„ - (5*) ' 



(26) 



Averaging Eq. ( p^ ) with the pro iahility distribution 
(EQ), one obtains the average LE, 



7 = In /a 



(27) 



Eqs. (|T]), (p^, and ( |2^ ) determine the two-point prob- 
ability distribution P2{zn, Zn+k) defined by Eq. (pjj[). 



B. Variance of the Lyapunov Exponent. General 
expression 

Eq. ( p^ for the variance of LE, cr^, can be presented 
in the following form: 



^ J2^Y1 (ln|2;„|ln|2;„+fc| 



(28) 



n=0 k=l 



+ i(ln^|z„|)-7^ 



The correlation function D{k) = (In |z„| In |z„+j,|) is in- 
dependent of the initial site n. With the use of the prob- 
ability distribution P2{zn, Zn+k) found in the previous 
subsection, it can be presented as 



oo 

„ , ImS f 



In |z| In 16(^)1 
(z-S) (z-S*) 



dz, 



(29) 



where ^fc(z) and 6 are defined by Eqs. (^4|) and (25), 
respectively. Interchanging the order of integration and 
summation in Eq. (p9) one can find for the variance 



_2 _ 2lmd f ln|z|ln|z5- 



{z^5) {z-5*) 



(30) 



- ^^(tt - ^) - |7ln {5^ - 1) -f 0(l/i2), 

where Lp is the phase of 5: 

6 = pexp(i(p), 



reduced to the interval [0, tt] and p denotes the absolute 
value of ^: 5 = \p\. The remaining integral in Eq. (|30| ) 
can be further simplified with the use of an appropri- 
ate contour in the complex z-plane. One, finally, arrives 
at the following expression for the variance in the case 
If < 7r/2, which corresponds to Uo > 0. 



dx tan ^ 



cosh(27a) — cos(2</?) 



sinh^(7a) 

sinh(27a) sin 1^9 
cosh(27a) cos — cos x 



(31) 

o(i/l2), 



Since our model is symmetric with respect to the trans- 
formation (/3 TT — (/3, the variance for (/3 > 7r/2 can be 
easily evaluated. 



IV. A NEW CRITERION FOR THE SINGLE 
PARAMETER SCALING 

The necessary (but not sufficient) condition for SPS to 
hold is that the localization length be greater than all 
microscopic scales in the system. Therefore, we should 
consider the general result, Eq. (|3l|), in the limit of large 
localization length, 7a <C 1. Our first goal is to develop 
an asymptotic form of the integral in Eq. ( |3l] ) in this 
limit. This is not a trivial task since the integral has 
a singularity at 7 = 0. The first term in Eq. ( ^ ) is 
also singular at this point, and one would anticipate two 
singularities to cancel out. The latter singularity has a 
logarithmic nature, 7 In 7, and we need, therefore, to ex- 
tract the similar logarithmic singularity from the integral 
in E q. (|3l| ). To this end, we first evaluate the integral in 
Eq. (|31|) by parts and present it in the following form: 



$(7, if) — TT tan 



-if P 



1 + C 

(C-cosip) 



^tan-if/^) (32) 



C — COS ip 



dx 



arccos {x + Q 



-i-C 



where parameters (3 and are defined as 
P = sinh(27a) simp, 



( — cosh(27a) cos p. 



(33) 



(34) 



Since we are interested in small values of 7a, we can ex- 
pand cos~^(a; -I- C) from Eq. (32) in the Taylor's series in 
x: 



cos ^(a:: + C)=cos -^(C) + ^a„a;" 



(35) 



n=l 



and carry out the term-wise integration. As the result 
we obtain 



5 



(36) 



tan 



1 



-1 ( C~cosip) 



In 



(C-C0S^)2 +/?' 



21 



(l + C)2 + /32 



+ /3[F(-(C-cos^))-F(-l-C)], 
where is defined by the foUowing power series 

r.n-1 



n=2 



1 



Comparing this expression with the original series ( |35| ) 
one can obtain for F{x) the following integral represen- 
tation 



1 1 / *N In Ixl 
F(.) = -eos-(C) + -^ 



The remaining integral in Eq. 
exactly. As a result we have 



j ^cos {C + x). 

(37) 

3^) can be calculated 



F{x) = - [cos- 1(C) - cos-i(C + x)] 
X ■' 

^ :ln|l-C'(C + a;) 



(38) 



+ V(l-C2)[l-(C + a;)2 



When one combines Eq. (gl|) with Eqs. (g7D and (|38|), 
the logarithmic singularity in $(7, f) nicely cancels out 
the singularity in the first term of Eq. (pi|). The ex- 
pression for the variance cr^ emerging in the leading in 7 
order takes the SPS form : 



27/L. 



(39) 



This is, to the best of our knowledge, the first truly mi- 
croscopic derivation of SPS with no ad hoc hypotheses. 
The main reward for this is the exact criterion for SPS, 
which follows from the conditions under which we have 
arrived at Eq. (^). First of all we assumed that 



■\/l — COS^ 



< 1. 



Since 



13 ~ 2-1 a sin If , C~l + 2(7a)^ 
the inequality (^) can be recast in the form 
7a ^ siniy9. 



(40) 



(41) 



(42) 



Another condition, which we have to impose in order to 
obtain Eq. (||) is 



(C — 1) cos(p/f3 ~ 07/ tani^ ^ 1. 



(43) 



Since sin ip < tan ip the first of the two inequalities is 
more restrictive, and the final condition for SPS takes 
the form 



hoc/Is > 1, 



(44) 



where koc = 7 ^ is the localization length, and a new 
length Is is defined as 



Is — a/ siniy9. 



(45) 



Evaluation of the integral in Eq. (31) in the limit 
K <ti 1 can be performed by means of a simple expan- 
sion of the integrand in power series in /3 and retaining 
only the linear in (3 term. The resulting expression for a 
can be presented as 



1 



21 



loc 



(46) 



Eq. ( ^ ) shows that in the regime considered here, is 
determined by the new length Ig rather than by lioc- It 
is important to emphasize that in this limit both lengths 
Is and koc can far exceed a, and Eq. (^), therefore, 
describes the violation of SPS while the system remains 
within a meaningful scaling regime. 

It should be noted, however, that Eq. ( p9| ) differs from 
Eq. (^ by the factor of 2. This discrepancy is due to the 
peculiar nature of the Cauchy distribution, whose mo- 
ments, starting from the second one, diverge. Because of 
this, none of the standard approaches, used to derive Eq. 
(^) within the random phase hypothesis, can be applied 
to the Lloyd model. In order to illustrate this point, let 
us consider, for example, an expression for cr^ obtained 
in Ref. M for AM: 



= —i Var 




In 1 



4-£;2 



(47) 



e e 
1 + " — f ns ;/ 



V4-£:2 




Brackets () designate here averaging over the random site 
energy e and the phase which is assumed to be statis- 
tically independent of e, and distributed uniformly. Var 
denotes the variance of the respective quantity. The stan- 
dard weak disorder expansion used in Ref. ^ implies the 
expansion of this expression in powers of the random vari- 
able e with consecutive averaging. The first term in Eq. 
( ^ ) then becomes of the order of ((e^) / [4 - E'^])'^ and 
is neglected, while the second term after averaging over 
the phase yields Eq. (^). In the case of the Cauchy dis- 
tribution for e, this approach cannot be applied because 
(e^) does not exit. In order to pass to the weak-scattering 
limit one has to average over e first, and only after that 
carry out expansion over the parameter P of the Cauchy 
distribution. Both terms in Eq. ( ^7| ) then become of the 



6 



same order of magnitude r/-\/4 — E^, and though the 
general proportionality between and LE is preserved, 
the numerical factor not equal to that of Eq. (||) or of Eq. 
(^). This result implies that the phase randomization 
hypothesis is not valid at all for the Lloyd model. What 
is more important, however, is the fact that although the 
phase randomization hypothesis fails, SPS still survives. 

For K <C 1, one can provide a clear physical interpre- 
tation for the length Ig. According to ThoulessJ^S the 
phase (p{E) is proportional to the integrated density of 
states G{E) = lp{E)/tt, and k ^ 1 corresponds to ei- 
ther G{E) < 1 or 1 - G{E) < 1. The length U then 
can be expressed in terms of the number of states in the 
energy intervals between E and the closest boundary of 
the spectrum Is = 1/{'kG{E)) [or 1^ = a/ [tt — airG)]. For 
AM these boundaries lie at ±oo, and for KPM they are 
the resonance boundaries of the bands, where i^{E) = 7rn 
with n an integer. The states in these regions arise due 
to rare realizations of the disorder, and can be associated 
with spatially localized and well-separated structural de- 
fects. The length Ig then can be interpreted as an average 
distance between such defects. In view of this interpre- 
tation of Is, the physical meaning of the transition be- 
tween two tjrnes of scaling regimes also becomes clear. 
Condition (|44| ) means that the localization length at the 
energy E exceeds the spatial separation between neigh- 
boring localization centers from the relevant part of the 
tail (between E and the nearest boundary of the spec- 
trum). Under this condition the localized states overlap 
and SPS is valid. 

To complete our discussion of the new scaling param- 
eter Is, let us compare it with the phase randomization 
length, Iph numerically studied by Stone, et al., Ref. |3^. 
Assuming that E is inside the conduction band and far 
from the band boundaries, we can approximate 1^ as 



0{Tl). 



(48) 



According to Eq. ([48|), Is decreases toward the center of 
the band C/q = 0, where it reaches its minimum value 
equal to a. At the same time, the phase randomization 
length was found in Ref. ^ to increase toward the cen- 
ter of the band E = 0, where it seemed to diverge. For 
E — 0, the phase distribution was found to be nonuni- 
form even for very long chains. The absence of the phase 
randomization in the center of the band was also found 
analytically in Ref. This comparison proves that Is is 
an independent new parameter responsible for the statis- 
tics of LE. Both numerical results of Ref. ^ and ana- 
lytical calculations of Ref. ^ show that a nonuniform 
distribution of phases can be consistent with SPS, pro- 
viding an additional argument against the condition for 
SPS based upon phase randomization. At the same time, 
our criterion, Eq. (^, correctly predicts validity of SPS 
in the band center as long as the localization length re- 
mains macroscopic. 



Is and the SPS criterion in generic models 

The peculiarities of the Lloyd model may cast doubts 
on the robustness of the new scale Is and the criterion 
Eq. (^). In order to show that this criterion is appli- 
cable beyond the Lloyd model, we carry out additional 
numerical simulations of the model studied previously in 
Ref. |3^. That model is of the Kronig-Penney kind, but 
unlike the model considered in the present paper, its po- 
tential is formed by rectangular barriers. The width of 
the barriers is assumed to be random with a uniform dis- 
tribution over a given interval. Both the potential and 
the statistics of the model used in numerical simulations 
is considerably different from the Lloyd model, e.g., all 
moments of the distribution function exist. 

It is instructive to rewrite expressions for in terms 
of a new dimensionless variable r = a^Llioc/2 as a func- 
tion of K. In terms of these variables both asymptotics 
of the variance, given by Eqs. ( ^9| ) and (p6|), can be 
presented in a form which contains no free parameters: 
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K < 1 
K > 1 



(49) 



Although we do not expect the concrete form of the func- 
tion t(k) to be universal, we do believe that the new 
crossover length Is retains its physical meaning in the 
general case, and that the crossover point is also univer- 
sally determined by k ~ 1. 

In order to generalize the crossover length Is for other 
models, we use the interpretation of the phase (p in terms 
of the integral density of states normalized in such a 
way that the phase would vary between and tt for any 
given band. The generalization is quite straightforward 
for models with a single band spectrum if in the absence 
of disorder the band has a finite width. The total number 
of states in such models is finite, and it can be used to 
normalize the phase. If the initial band of the system is 
infinitely broad, e.g., for the Schrodinger equation with 
a random white-noise potential, one has to introduce a 
cut-off frequency for the spectrum in order to normalize 
the phase. The crossover parameter k for the later model 
can be obtained from an expression relating LE and the 
integral density of states found in Ref. ^: 



IX exp 
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Ex 



dx, (50) 



where D determines the strength of the (5-correlated po- 
tential. The asymptote of k for large negative E is 



exp 



A\E 



,3/2' 



3D 



« 1, 



and we conclude that these states do not obey SPS. Tran- 
sition to SPS behavior again occurs at the initial band 
boundary E = 0, where k « 1.1 and does not depend 
upon parameters of the model. 
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For systems with multiple (in the absence of disorder) 
bands separated by band-gaps, one has to consider sepa- 
rately two different situations. If disorder-modified bands 
still do not overlap and a genuine gap between the bands 
persist, the situation is equivalent to the single-band case. 
The phase can be defined for each separate band and nor- 
malized by the number of states in the band. The results 
of the numerical simulations of this particular situation 
are shown in Fig. 1 along with t(k) obtained from our 
analytical Eq. Is was calculated numerically with 

the phase defined using the integral density of states nor- 
malized in such a way that the phase changes from zero 
to TT when energy sweeps over a band from one fluctua- 
tion boundary to the other. One can see from this figure 
that the crossover between different asymptotes for both 
numerical and analytical calculations occurs in the same 
region. This proves the universal significance of our cri- 
terion (44) for SPS and justifies the suggested general- 
ization of the crossover length Ig. 



1.20 



1.00 - 



0.80 



0.60 



0.40 



0.20 



0.00 




region where the tails from different bands overlap). Let 
gn (E) be the differential density of states related to the 
n-th band. Then integral Ntot = 9niE)dE gives the 
total number of states originating from the band. One 
can define a phase 



ipn{E) = TT / 



gn{E)dE/Ntot, 



with E obeying the following inequality i?,"7n < E < 
■^mim where E^^^ corresponds to the minimum of the 
actual differential density of states within the gap be- 
tween the n-th and (n -|- l)-th bands. The phase defined 
according to this procedure does not assume unphysi- 
cal values of tt at the points where there are no spectrum 
boundaries, and we suggest that the parameter Ig defined 
through this phase according to 



L ^ = sin 



Nt, 



gn{E)dE 



rnin - min'> 



(51) 



can be used in order to formulate the criterion Eq. ( |44| ) . 
The suggested definition of Is can be practically used for 
analytical estimates of the transition between different 
statistics, using, for example, a tight-binding approach to 
a multi-band problem, where interactions only between 
adjacent bands arc taken into account. However, more 
detailed discussion of this issue requires a separate paper. 



V. PROPERTIES OF THE TRANSITION 
REGION BETWEEN SPS AND NON-SPS STATES 



FIG. 1. The function t{k) obtained from the analytical so- 
lution of the Lloyd model Eq. ( ^ ) (solid line) and the numer- 
ical simulations of KPM with rectangular barriers of random 
widths (open circles). Note that the SPS equations for these 
models differ by a factor of 2 [Eqs. (^) and (^9|)]; therefore 
we rescaled numerical data by this factor. 

A new situation arises, however, when fluctuation 
states from adjacent bands overlap, and the spectrum 
does not have boundaries in former band-gaps. On one 
hand, it is clear that if the number of states in the for- 
mer gap is small, SPS should not be expected on the basis 
of the general qualitative interpretation of the criterion 
(^). On the other hand, since there are no exact bound- 
aries of the spectrum inside former gaps, one cannot de- 
fine a phase suitable for determination of Is- Though 
this situation requires special consideration, we can offer 
a conjecture that can be used to meaningfully define Is in 
this case. Consider one of the original bands between two 
adjacent band-gaps. In the presence of disorder there ap- 
pear tails of the density of states within the gaps. When 
disorder is small, one can always distinguish between gap 
states originating from different bands (except for a small 



In this section we discuss properties of the transition 
region between SPS and non-SPS states. In spite of the 
mentioned peculiarities of the Lloyd model, our calcu- 
lations provide a sound qualitative explanation for nu- 
merical results of Ref. |3^, confirming once again that we 
correctly describe the qualitative nature of the transition 
between SPS and non-SPS statistics. In the model con- 
sidered, the phase, tp, and LEL 7, can be conveniently 
presented in the following formQ 



sm if 



(52) 



(53) 



where s is given as 
1 



(4 + r: 



u 



(4- 



^ u 



(54) 



and we assume again that 7^1. The relation between 
7 and if is determined by the parameter s, which in its 
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turn depends upon the energy E. Let us recall that the 
energy enters into our equations through parameters Uq 
and Tjj defined in Eqs. (||) and (^). For energies within 
conduction bands, LE is of the order of Tu though Ig is 
of the order of one [sec Eq. (H)]. Thus, SPS holds as 
long as disorder ispsmall, Tij <C 1, in accordance with the 
previous results MeiyEl From Eqs. (|5|), (H), and (H) 
it is clear that the relation between lioc and Is changes 
with the energy approaching an initial spectral bound- 
ary. In the limit of small disorder koc = h exactly at the 
boundary. Therefore, one should expect the strongest vi- 
olation of SPS for states which arise due to disorder in the 
originally forbidden regions. For AM this corresponds to 
energies \E\ > 2, and for KPM these are energies from 
band-gaps of the original spectrum. For energies liying 
far away from the boundaries one can obtain the follow- 
ing approximate expressions for LE and the phase ip. 



Tr 



sm if 



7 ' 



Un 



«1, 



(55) 
(56) 



It is evident that in this case hoc <C Is, and the variance 
behaves according to Eq. (^6|). The states disobeying 
SPS, however, are more important for KPM than for AM. 
The reason for this is that LE in AM becomes of the order 
of one not very far from the boundary, moving the system 
out of any scaling regime. In KPM 7 can remain small 
enough throughout entire band-gaps for sufficiently high 
energies k ^ Vq, and the violation of SPS in this case oc- 
curs when the localization length is still of a macroscopic 
scale. 

Eqs. ( |39| ) and ( ^ ) explain a non-monotonic behav- 
ior of cr{E) observed numerically in Ref. When the 
energy moves towards a band edge, LE grows and a 
grows along with 7. When, however, 7 becomes equal to 
Ij^ the variance, starts decreasing towards the value 
~ ttIj^/N. The maximum of cr, therefore, corresponds 
to the energy where 7 ~ Ij^, i.e. the boundary of the 
original spectrum. 

We can now also estimate the width of the transition 
region between SPS and non-SPS states, which was found 
in numerical simulations to be surprisingly small. The 
transition between the two groups of states occurs when 
(j{E) passes through its maximum, and the width of the 
transition region is related to the sharpness of the max- 
imum. In view of the preceding discussion, the latter is 
determined by the region of energies over which 
changes its behavior. The extent of this region can be 
estimated from the condition 



"0 1 2C^or[/ 

In AM it leads to 6E ^T, in KPM one has 

r 



6k . 



Vo 6n ' 



(57) 



(58) 



where A„ is the width of n-th band-gap, fc„ represents 
the non-resonant boundary of the n-th band, and the 
parameter Sn is defined as 



A 1 ^ 



(59) 



In both cases the width of the transition region is de- 
termined by the degree of disorder in the system, and is 
small when disorder is small. In the Kronig-Penney sit- 
uation, however, Eq. (|5^) indicates a special sharpness 
of the transition in the case of high-energy bands, when 
the parameter A„ is also small. 

When disorder increases,pAM and KPM behave differ- 
ently. Monte Carlo resultsE3 show that in periodic sys- 
tems, an increase of disorder leads to a restoration of SPS 
for almost the entire spectrum of the system. We are now 
able to explain this behavior and to provide an estimate 
for the critical disorder. It is clear that the parameter 
Ij^ reaches a minimum at the energy in the center of a 
band gap. This minimum value can be estimated from 
Eq. (M) as 



r 

kn 



(60) 



where fc„ represents the non-resonant boundary of the n- 
th band, and the parameter Sn is defined by Eq. ( |59| ) . At 
the same the energy where is minimal, LE assumes 
its maximum value: 



^^0 



(61) 



Cmin increases with disorder, while 7inax does not change, 
and at F ~ Vo/Sn two length scales are of the same or- 
der, Zsmin — 7max- this instant for the states outside 
of the immediate vicinity of the center of the band-gap, 
Is ^ 7~^ and SPS is restored. Thus, we can identify 
F"^ — Vo/Sn as a critical disorder for the n-th band-gap. 
For the states right in the center of the gaps, however, 
h ^ 7~^ no matter what disorder is, and these states do 
not obey SPS. Therefore, the complete restoration of SPS 
for the entire band-gap does not occur in this model, but 
the width of the non-SPS region decreases with increase 
of disorder, as one can see from Fig. 2. 
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nn, both the locahzation length koc and the crossover 
length Ig diverges, while the variance of LE vanishes. Al- 
though the resonances are not stable with respect to a vi- 
olation of the periodic arrangements of the (5-potentials, 
they occur in some other models as welL_. We already 
mentioned models with correlated disordertZl and random 
superlattices.EaS The latter has an experimental signifi- 
cance with applications to propagation of classical waves. 
Therefore, it is interesting to consider the behavior of the 
critical parameter k in the vicinity of the resonances. Al- 
though both lioc and 1^ diverge at the resonances, their 
ratio K remains finite and takes on the following values: 



31.5 32.0 32.5 33.0 33.5 34.0 34.5 



To 

n 

To 



1 + ^, ka<.n 
To 



— . ka > TTn 
To 



(62) 
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One can see from Eq. ( |62| ) that k experiences a disconti- 
nuity at resonance points: its value decreases by 2Vo/To 
once a point is crossed. In the case of small disorder, 
when Vq/Tq ^ 1, this is a dramatic jump, such that 
K ^ 1 at the band side of the resonant boundary and 
K <C 1 at the gap side. It is obvious, therefore that scal- 
ing properties of the system also change discontinuously 
at the resonance from SPS behavior at the band side to 
the scaling with two parameters at the gap side. 



FIG. 2. The variance of LE (solid line) and LE itself 
(dashed line) for KPM are posted as functions of the energy 
parameter, k, for different degrees of disorder. The plot ex- 
tends over the region lOvr < k < IItt, so one gap and one 
conductance band are represented. The presence of the max- 
imum in the vicinity of the band boundary {k ~ lO.SSvr) sig- 
nals the violation of SPS for states from a band-gap. The 
absence of a similar maximum at the second boundary is due 
to its resonance nature. The degrees of the disorder, F, are 
3 (a), 7 (b), and 15 (c). It is clearly seen that with the in- 
crease of the disorder, the number of states disobeying SPS 
decreases. 

In the case of a strong potential Vq and in the high 
energy, <^Vo<^ kn, one has Sn — 1, and the critical 
disorder is just equal to Vq. When a potential is weak, 
Vqcl <^ 1, there exists a medium-energy regime, when 
i/Vb/a ^ kn ^ Vq. In this case 

\/^)a^ 1, 
/'I - 

'smin ^ ^2 Y a ' 

It is interesting to note that in this case F"^ increases 
with an increase of kn, even though the widths of the 
gaps decrease. 

The calculations presented above referred to the non- 
resonant band boundaries. At the resonant points, ka = 



VI. CONCLUSION 

In this paper we studied statistical properties of the 
Lyapunov exponent in the one-dimensional Anderson 
model with the Cauchy distribution of site energies. The 
model can also be interpreted as the Kronig-Penney 
model with periodically positioned (5-potentials with ran- 
dom strengths. The main objective of the study was to 
find an exact solution for the thermodynamical limit of 
the variance of LE and to establish an exact criterion 
for the existence of single parameter scaling. It is im- 
portant to emphasize that in contrast with all previous 
calculations of the variance, we did not use the phase 
randomization hypothesis. This allowed us to reject the 
generally accepted assumption that it is the length over 
which the phase of reflection and transmission coefficients 
becomes uniformly distributed that sets the condition for 
the existence of SPS. 

We found a new length scale, Ig, which is responsi- 
ble for the scaling properties of the conductivity in the 
system: SPS exists as long as the localization length 
koc exceeds Is- The length however, differs from the 
phase randomization length, and presents, therefore, a 
new significant scaling parameter. The parameter 1^ is 
microscopic for states close to the center of the original 
conduction bands of the system and does not impose, 
therefore, any additional restrictions for the existence 
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of SPS excepting the regular requirement for the local- 
ization length to be of a macroscopic dimension. How- 
ever, for the states at the edge of the bands, Ig grows 
to a macroscopic size and the condition lioc — Is actu- 
ally establishes a boundary between the states with and 
without the SPS statistics. As soon as Is becomes much 
larger than all microscopic lengths, this scale becomes 
significant. In this limit it can be expressed in terms of 
the number of states Ndef{E) which arise at the tails of 
the initial bands due to rare fluctuation configurations: 
Ij^ = Ndef{E). It then can be given a natural phys- 
ical interpretation as an average distance between such 
defects. 

The change of the scaling behavior occurs when the 
energy crosses over a boundary of a former gap. In the 
case of regular boundaries, the change occurs gradually 
with the critical parameter k being of the order of unity 
right at the boundary. The Kronig-Penney version of 
our model, besides regular boundaries has so called res- 
onant boundaries, where both LE and l^^ vanish. We 
found that at the resonance boundaries the parameter k 
undergoes a sudden jump from very large values k ^ 1 
at the band side of the boundary to very small values 
K <C 1 at the gap side. This means that the change of 
the scaling behavior at the resonant energies also occurs 
discontinuously: the system obeys SPS when the bound- 
ary is approached from the conduction band, and then 
demonstrates two parameter scaling if the boundary is 
approached from the gap. 

We carried out numerical simulations of the Kronig- 
Penney-like model with a different configuration of the 
potential and different statistics. The comparison be- 
tween numerical and analytical results clearly indicates 
that significance of the length scale Is defined in terms of 
the integral density of states persists beyond the Lloyd 
model, and that the new criterion for SPS established in 
the present paper has a universal nature. 
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